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This article examines the time-dependent Hartree-Fock (TDHF) approximation of 
£f~^ , single-particle dynamics in systems of interacting fermions. We find the TDHF approx- 

imation to be accurate when there are sufficiently many particles and the initial many- 
particle state is a Slater determinant, or any Gibbs equilibrium state for noninteracting 
fermions. Assuming a bounded two-particle interaction, we obtain a bound on the error 
of the TDHF approximation, valid for short times. We further show that the error of the 
the TDHF approximation vanishes at all times in the mean field limit. 



1 Introduction 



Dirac invented the time-dependent Hartree-Fock equation in 1930. The time-dependent 
Hartree-Fock (TDHF) equation is a nonlinear Schrodinger equation designed to approximate 
the evolution of the single-electron state in an n-electron system. Dirac noted that the TDHF 
equation, originally written as a system of n coupled Schrodinger equations for occupied orbitals, 
may be written as a Liouville-von Neumann equation for the single-particle reduced density 
operator. We study the TDHF equation in this form, availing ourselves of trace norm techniques 
inspired by ^2] to estimate the error of the TDHF approximation. 

We will consider Hamiltonian dynamics of fermions interacting through a two body potential 
V . The energy operator for a solitary particle will be denoted by L, the energy of interaction 
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of a single pair of particles will be denoted by V, and the total energy operator for a system of 
particles will be the sum of all single-particle energies and all pair energies. In this article we 
only consider bounded interaction potentials V; the case where V represents Coulomb repulsion 
between electrons will be treated in a forthcoming paper [3] . Although the number of particles 
does not change under the dynamics just described, we prefer to formulate the dynamics on a 
fermion Fock space, so that we may consider initial states of indeterminate particle number. We 
are particularly interested in initial states which are Gibbs equilibrium states (grand-canonical 
ensembles), for we are going to show that such initial states enhance the accuracy of the TDHF 
approximation. 

Let L be a self-adjoint operator on HI, and let V be a bounded Hermitian operator on H®H 
that commutes with the transposition operator U defined by U(x £g> y) = y <g> x. We are going 
to discuss the dynamics on the fermion Fock space Fh whose Hamiltonian H may be written 
in second quantized form as 

H = )]{j\L\i)ajai + (kl\V\ij) a\a\ajai (1) 

i,j i,j,k,l 

We will analyze the solutions of the Liouville-von Neumann equation 

ihjD{t) = [H,D(t)\ 

D(0) = A), (2) 

which is the evolution equation for the density operator in the Schrodinger picture of quantum 
dynamics. 

We will see that (j2J leads to the following equation for the single-particle density operator 

M(0) = Mpo)- (3) 

(The subscript : i denotes a partial trace; see definition (fTT| below.) This equation for Mi(t) 
is not "closed," for its right hand side involves the two-particle density operator M^if)- The 
TDHF approximation to M\(t) is the solution of the initial value problem 

ihjF{t) = [L,F(t)\ + [V, F(t) 82 2A 2 ] :1 

F(0) = M(5(0)) (4) 

where A 2 is the orthogonal projector of EI ® EI onto the subspace of antisymmetric vectors. The 
existence and uniqueness of solutions of (j3J) were established in jB] for the case where V is a 
bounded operator, and in Ej for the case where V is a Coulombic interaction. 
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The TDHF equation (J1J is obtained by closing the single-particle equation J3J) with the 
Ansatz 

■A/2 = (Mi ®Nx)2A 2 (5) 

at all times. The relation (jSJ) holds for pure states corresponding to Slater determinants, and 
also for Gibbs densities. However, even supposing that A/2(0) satisfies ©, the interaction V 
is likely to introduce "correlations" in J\f 2 (t), that is, departures from (j5J), and ignoring those 
correlations in the TDHF equation requires some justification. 

We are going to prove that the absence of correlations is self-perpetuating in the mean field 
limit: if A/^O) satisfies © then Mz{t) asymptotically satisfies (0) as the number of particles 
N tends to infinity and the interaction strength is scaled as 1/N. This is the content of 
Theorem 15.21 To prove this theorem we first bound the error of the TDHF approximation 
in terms of the average particle number and interaction strength. This bound is presented 
in Theorem 15.11 without reference to the mean field scaling, and it applies to any system of 
fermions whose Hamiltonian has the form (Q). Unfortunately, the bound is valid only at short 
times, i.e., up to a time inversely proportional to || V]| and the average particle number. 

Let us advert to some shortcomings of our results. Firstly, we consider only bounded two- 
particle potentials V. Fortunately, the challenges presented by the Coulomb potential can 
be overcome j3[ E] and the mean field limit of Theorem 15.21 also holds for certain electronic 
systems. Secondly, the explicit error bound of Theorem 15.11 is valid only at short times, which 
are too short to be of interest in molecular-electronic problems, even if the Coulomb interaction 
between electrons is truncated at the Bohr radius. However, the model (0) does not only apply 
to molecular-electronic problems, and we hope that Theorem 15.11 may find other applications, 
perhaps to certain models of interactions between nucleons j3]. 

We first published our derivation of the TDHF equation in the mean field limit in [2|. There 
we assumed that the initial states are Slater determinants. The main theme of this article is 
that the initial states need not be Slater determinants; the TDHF approximation should work 
equally well (or badly) for all initial many-fermion states of Gibbs type. Also, in this article 
we are not only interested in the mean field limit, and we offer the error bound of Theorem 15.11 
for the unsealed problem. 

The next section gives the background on fermion Fock space, trace class operators, reduced 
density operators, and Gibbs equilibrium states. Section |3] carefully defines the many-particle 
dynamics generated by (0) and derives the single-particle equation Q for A/i(t). Section |U 
introduces the TDHF approximation of A/i(i) and its higher-order analogs. Section presents 
Theorems 15.11 and 15.21 Finally, the appendix contains the proofs of several Propositions. 
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2 Definitions and notation 



Consider a quantum particle whose Hilbert space is HI, i.e., a particle which, in isolation, would 
constitute a system whose (pure) quantum states are represented by the rank-one orthogonal 
projectors on some Hilbert space HI. The set of quantum states available to a system of n 
particles of this kind is determined by their "statistics," i.e., whether the particles are fermions, 
bosons, or distinguishable. If the particles are fermions, the pure states of a system of n of 
them are represented by the rank-one projectors onto vectors in the antisymmetric subspace 
H[( n ) of the tensor power space H® n . 

The Hilbert space M® n is the closed span of the simple tensors 

xi © x 2 © • • • © x n x\, x 2 , ■ ■ ■ , x n G H 

with the inner product 

(yi ® V2 © ■ ■ ■ © y n , xi ® x 2 ® • • • ® x n ) = (yi, x{) (y 2 , x 2 ) ■ ■ ■ (y n , x n ) . 

Let n n denote the group of permutations of {1,2,..., n}. For each n e U n , a unitary operator 
U n on H® n may be defined by extending 

U^Xi ® X 2 ® ■ ■ ■ ® X n ) = 2^-1(1) ® ^-1(2) <8> • • • ® V!(n) 

to all of H ". A vector ip e H® n is antisymmetric if U n (ijj) = sgn(7r)^ for all 7r e n n . 
The antisymmetric vectors in HP" 1 form a closed subspace which will be denoted rf™). The 
orthogonal projector with range is 

A « = A S sgn(vr)^. (6) 
n! ^— ' 

7ren n 

If Xi, X2, . . . , x n is an orthonormal system in a Hilbert space HI, then the vector 

Vn! A n (xi (g x 2 <g> • ■ ■ ® x n ), (7) 

is a unit vector in M^ n \ A vector of the form (j2J) is called a Slater determinant. If {e a } is 
a complete orthonormal system in HI, then a basis of HI (n ) may be chosen from the set of all 
Slater determinants of the form (JJJ) where {xi, . . . ,x n } a subset of the single particle basis of 
cardinality n. 

If the number of particles in the system is not fixed, the appropriate Hilbert space is the 
direct sum of the n-particle spaces M^ n \ This is the (fermion) Fock space 

F e = H (0) © H (1) © H (2) © H (3) © • • • . (8) 
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The possibility of a zero-particle state is accommodated by m C, a one-dimensional space 
spanned by the vacuum vector Q. The number operator N on Fe is the self-adjoint operator 
whose restriction to H^ n ) equals multiplication by n. Annihilation and creation operators a x 
and a\, are explicitly represented on Fe as discussed in |S]. The statistical state of a system of 
fermions with single-particle space H determines a bounded positive continuous functional uj 
on the bounded operators on Fe with = 1. We are only going to consider normal states, 
i.e., those uj that satisfy 

oj(A) = Tr(D u A) 

for some nonnegative trace class operator D u of unit trace, which may be called the (statistical) 
density operator. 

Density operators on any Hilbert space IK are trace class operators in particular. Let T(K) 
denote the real Banach space of Hermitian trace class operators T with the norm ||T||i = 
Tr(|T|). We often use the inequality that, for bounded operators B on K, 

\Tr(TB)\ < ||T||i||B||. (9) 

A linear functional defined on the space of compact Hermitian operators by 

$t(K) = Tt(TK) (10) 

is continuous by Q, and indeed the space T(K) is isomorphic to the Banach dual of the space 
of compact Hermitian operators on K, via the isomorphism T i — > 



2.1 Reduced density operators 

We now restrict our attention to density operators D on Fe that commute with the number 
operator and such that N m D is trace class for all m G N. For such D we will define the m th 
order reduced number density operators j\f m (D) and explain their physical significance. 

We begin by defining the reduction of an n-particle statistical density operator, i.e., a positive 
operator of unit trace D on M^ n \ to an m-particle density operator denoted D. m . If A is an 
operator on H.( n \ let A denote the extension of A to H®" defined by first projecting onto HI*™) 
and then applying A, i.e., 

A = A © Oj|®riQ]jj(n) . 

For each < m < n there exists a positive contraction from T(M^) onto T(H^ m ^) known as 
the partial trace. This partial trace map, which we denote by T i— > T. m , is defined implicitly 
through the duality (jl(Jj) by the requirement that 

Tr((if(g)/® ...® I)T) = Tr(KT.. m ) 
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for all compact Hermitian operators K on HI (m < ) . It follows that, for any orthonormal basis O 

of e 0n - m , 

(£,r :m V> = ^((£©o)®0, r((^©o)®0)> . (n) 

In case D is a density operator on EU n ), and if m < n, the operator P :m is known as the 
m-particle reduced density operator [TT]. It is used to determine the expected values of m- 
particle observables when the n-particle system is in the statistical state D. The m-particle 
reduced density operators obtained from P^ will be denoted (P^,) . If ip is an n-particle Slater 
determinant, then the reduced density operators (P^). m satisfy 



1 n 

(p,) :1 = -£p, 



n z — ' 

i=i 

(^) :m = (13) 

Now we can define the reduced number densities jV m , which serve to describe the m-particle 
correlations in a system of many particles. Let D be a density operator that commutes with 
N. Then 

CO 

D = 0A, (14) 
n=0 

where each P n is a nonnegative trace class operator on H^ n ). Assume that 

CO 

J2n m Tr{D n ) < oo. (15) 

n=0 

For such D, define the m th order reduced number density 



CO . 



M m {D) = Y j , — -D n:m . (16) 

z — ' [n — my. 

n=m v ' 

The operators J\f m (D) are called reduced density operators in Section 6.3.3 of JH], but we prefer 
to call them reduced number density operators because the trace of M\ is the average particle 
number. They are trace class by (fT3j) . An mparticle density operator D n extends to a density 
operator 

D n = 0©---©0©P n ©()©■■• (17) 

on all of Fock space. It is clear that Af m (D n ) equals Af m (D n ) = n(n — 1) • • • (n — m + l)D n:m if 
m < n but it equals the zero operator if m > n. 

Formula (|19jl below may clarify the sense in which M\ determines expected values of single- 
particle observables. Following [H], we define the second quantization dT(H) of a self-adjoint 
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operator if on HI to be the closure of the essentially self-adjoint operator (BdT n (H) on F H , 
where dT n (H) denotes the restriction to of the operator 

n j—l times n—j times 

dT n {H) = 22 I® I® •-•®J (18) 

on H® n . The operator dY(H) is called the second quantization of H; we might also call it a 
single-particle observable. Suppose that u is a normal state on Fe whose density D u satisfies 
flUJ) and (HSJ). From (JTHJ) it then follows that 

u(dT(B)) = Tr(N"x(D u )B) (19) 

for any bounded Hermitian operator B on EL For example, the number operator iV on Fe is 
the second quantization of the identity operator, i.e., N = dT(I). From (JT§j) we see that the 
trace of M\{DJ) is the average number of particles in the system when it is in the state to. 
Another example: if x is a unit vector in H, the operator a x a x on Fe is the second quantization 
dT(P x ) of the rank-one projector P x , and (jTUj) tells us that oj{a\,a x ) = Tr(J\fi(D u )P x ). 

Theorems 15 . II and 15 . 21 relv on the following fact, which is proved in the appendix: 

Proposition 2.1 If D is a density operator on Fe that commutes with N and such that 
Tr{ND) < oo, then the operator norm of Ai(-D) is not greater than 1. 

2.2 Gibbs equilibrium states 

A Gibbs equilibrium state is that of a system of noninteracting fermions in thermal and chemical 
equilibrium with its environment, but only in case e~^ H is trace class, where H is the single- 
particle Hamiltonian and 1/(3 is Boltzmann's constant times the temperature. In this case 
the density operator is proportional to exp(— ft dT(H — ///)), where \i G K is the chemical 
potential. This density operator is diagonalizable with respect to an occupation number basis 
of Fock space, and occupation numbers are stochastically independent. If the single-particle 
Hamiltonian has discrete eigenvalues ei < e 2 < • ■ ■ , the probability of occupation of the j th 
level equals 1/(1 + e^~^). 

When D is a Gibbs equilibrium density then J\f m (D) = Afi(D)® m m\A m . The same is true 
when D is a pure state corresponding to a Slater determinant. Indeed, this relation characterizes 
the "gauge-invariant quasifree states of the CAR algebra" that have a trace class single-particle 
operator, and both Slater densities and Gibbs equilibrium densities are of this type. 

We describe these states in probabilistic language. Let {4>j}j^j be a basis of EL One basis 
of Fe consists of the vacuum vector Q and all Slater determinants made of vectors from {4>j}- 
This basis is indexed by J 7 , the set of finite subsets of J including the empty set. A nonempty 
subset s = {ji, . . . ,j n } of J corresponds to a Slater determinant \I/(s) formed from the vectors 
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(f)j 1} . . . , <j)j n . (There are in fact two such Slater determinants, opposite in sign, and we choose 
one of them.) The empty subset of J corresponds to the vacuum vector Q, i.e., P(0) = fi. Let 
P be a probability measure on T with the cr-field of all its subsets. That is, let P : T — ► [0, 1] 
be equal to except on a countable subset of J 7 , on which subset J^P(s) = 1. Define the 
random variables Nj on T by 

, r / N r 1 if j 6s 
N ^ = { if U b 

and define N — ^2Nj. Suppose that the Nj are independent with respect to P, and define 

P(JV i = 1) = p(j). (20) 

Then 

n 

P(s) = P((% = i) n • • • n (N jn = l)) = JJpta) (21) 

83j'l,...,jn i=l 

when ji, . . . ,j n are distinct. Note that E(iV) = p(j) < oo since P(iV < oo) = 1. 
Given a probability measure P on J 7 , we define a density operator on Fh by 

G[P] = ^P(s)P s , (22) 

where P s denotes the orthogonal projector onto the span of ^(s). 

Proposition 2.2 Let P be as in and WIS) . Let G denote the density operator G[F] of Wty . 
Then for all n GN 

N n {G) = Ni{G)® n n\A n . (23) 

The preceding proposition is proved in the appendix. Proposition 12.21 and the following propo- 
sition (also proved in the appendix) together imply Proposition 12.41 below: 

Proposition 2.3 If T is a Hermitian trace class operator then 

\\T® n n\A n \\ l < \\T\\[. (24) 

Proposition 2.4 Let G be as in 

Proposition^ Then (^(G)^ < ||A/"i(G)||". 
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3 Definition of the dynamics on Fock space 

We now define the dynamics (0) more carefully, and derive the reduced dynamics (jHJ) from (J2J). 

Let if be a self-adjoint operator on H, and let V be a bounded Hermitian operator on EI <S> HI 
that commutes with the transposition operator U(x®y) = y <8> x. For 1 < j < n, let Lj denote 
the operator 

j— 1 times n—j times 

I <g> • • • ® I ®L® I <g) • • • ® I 

on H® n (the value oi n > j is not explicit in the notation Lj but it will always be clear from 
context). For 1 < i < j < n, let L/^ denote the permutation operator on H®" that transposes 
the i th and factors of any simple tensor x\ <g) ■ ■ ■ ® x n , and let 

Vii = u m u (2j) {V ® I® n ' 2 ) u (2j) u m 

(again, the domain M® n of Vij will always be clear from context). For each n, define the 
operators 

n 
3=1 

H (n) = L (n) + Vy 
l<i<7<n 

on (these operators are defined on all of M® n but we are only considering their restrictions to 
the invariant subspace H^). The Hamiltonian operator H, which we had formally represented 
above by (^J), is the direct sum H = ® defined on the domain 

V(H) = 1 2 = ®x n G Fe : ^ \\H^x n \\ 2 < oo}. 

n 

This operator is closed and self-adjoint (see Section 6.3.1 of 8j), and — jM is the generator of 
the strongly continuous group 

oo 
n=l 

of unitary operators on Fh, where = exp ( — jrH^). 

It is convenient to have some notation for the free part of the dynamics, and we will subse- 
quently use 

. oo 

U} n) = exp(-|i») and U t = 0f/ t W . 

n=l 
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The Liouville-von Neumann dynamics corresponding to (J2J) are given by the group 

W t (D) = W t DW- t (25) 

of isometries of T(Fe). (See Proposition 3.4 of [0] for a proof that groups of isometries defined in 
this way are strongly continuous.) Define the subspaces T n C T(H) consisting of all n-particle 
trace class operators: 

% = {fas in (HID : r n 6 T(H^)} . 
These subspaces are invariant under W*, and the restriction of Wt to 7^ is 

wi n \T) = wl n) TW%\ 
The generator of this group of isometries is 

l<i<j<n 

where is the the generator of the group 

l4 n \T) = lll n) TU!$, (27) 

n 

which may be written formally as — \ ^2 [Lj , ■ }. (See jH] and references therein for a proof.) 
Since (J2fij) is a bounded perturbation of C^ n \ it follows that W t (T n ) equals T n (t), where 



T n (t) = ut ] T n (0) - % - \ U£\ lV^T n (s)]ds 

l<i<i<n 

when T n E T n . Taking the m th partial trace of both sides of (|2S|) yields 

i /"' 



(28) 



l<i<j<m 



— (n — m)-j: 



j'=i 



because of the symmetry properties of T and V. Multiplying both sides of the last equation by 
n\/{n — m)\ we obtain 

-r.Tn-.mit) = U[ ^ rrT n:m (0) — — / [v^-, rr^n:m(s) 

[n — my. (n — my n J n , jr—f^ L m — ml! 



l<«<j<m 



- ^TyT„ :m+1 ( S )j d a . (29) 

" ~^ L (n — m — 1J! 
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Now let D = ®D n be a density operator on F H . Then W t D = Q)D n (t), where D n {t) is the 
solution of (J28j) with initial condition D n (0) = D n . If D satisfies moment condition (|15|) then 

oo . 

E71 ' 
7 '-^D{t) nm 
(n — my. 

n=m v ' 

is trace class. We abbreviate Af m (D(t)) by Af m (t). Summing the right-hand sides of (|2T?j) with 
D in place of T yields 

pt 

N m {t) = ui m W m (o) - \ \ utl Wa^M] 

l<i<j<m 

t m 

«S5 E [^,- + i,AT m+1 ( s )] :m rf s , (30) 

for (|15p permits the interchange of the sum and the integral. This is equation © in integral 
form when m — 1. To summarize: 

Proposition 3.1 Let W t and Wt 6e as defined in \27\) and 

Suppose that D is a density operator on Fh of the form D = ®D n such that moment condition 
/ fT3)) holds for some m G N. LetJ\f m (t) denote J\f m (Wt(D)) . ThenJ\f m (t) satisfies equation IRJIty. 




4 The TDHF hierarchy 

The existence and uniqueness of mild solutions of the TDHF equation (£Q) is established in j^J. 
There it is shown that the integral equation 

F(t) = U^F(0) - y^i-l[V, F(sr 2 2A 2 ] :1 ds (31) 

has a unique solution F(t) for any Hermitian trace class operator -F(O). Define T\{t) = F(t) 
and, for m > 1, define 

F m (t) = F(tf m m\A m . (32) 

We proceed to derive equations for the J r m (t) from (|3*T|. 
First, set G(t) = U^F(t), so that 

[V, F( S )^ 2 2A 2 ] :1 d S . 
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Now apply the product rule (in integral form) to G(t)® m : 

m rt 



G(tf 



3=1 
-t 



F{0f m - -J2 I Gis)®*- 1 ®US [V, F(s)® 2 2A 2 ] :1 ®G(s)® n - j ds 



F(0f m - - I U 



(m) 



m 



ds. 



3=1 



Apply to both sides of the preceding equation to obtain 
F(t)® m = ul m) F{Q)® m - 



ft m 

/ Z42£[W> F{sT m +\l-U {hm+1) ) 

J ° 3=1 



ds. 



Multiply both sides of the last equation by m\A m on the left, noting that lis {X)A r< 
Us m \XA m ) and that A m commutes with ^2 Vj, 



" j,m+l • 



T.M) = ut m) F m (0) - I [«SE[W F{sr m+ \l-U (hm+l) )m\A r 



3=1 



ds. 



Since (m + l)\A m+ i = (i — U { 



(l,m+l) ' ' ' 



U( mjm +i))m\A m , the last equation may be rewritten 



F m {t) = U^^) - - / ZYS^[^ +1 , F m+ i(s)]. m ds 

n J ° 3=1 



l<j7^A:<m 



M 

t-s 



V,- m+1 , F(sf m+1 f/ (fc , m+1) m!A 



(33) 



(is. 



We call these equations for the T m {t) the TDHF hierarchy. 
The trace norm of the last term in is bounded by 

2 /"* 

m(m-l)- J {V m -l >m +lU(rn,m+l){Fm{s) ® F(s))} , 

It can be verified that 

{V m ^ m+1 U {mim+1) (^ m (s) ® F(s))) :m = (J ™- 1 ® F(s))y m _i, m ^ m (s), 
whence the trace norm of the last term in ()33|) does not exceed 



m{m-l)lj\\V\\ \\F{s)\\ \\F n {s)\\ x d8. 



(34) 
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Now ||.F TO (s)||i < \\F(s)\\™ by Proposition O since T m = F® m m\A m . Furthermore, \\F(s)\\x = 
\\F{0)\\! and ||F(s)|| = ||F(0)|| for all s > by Proposition 4.3 of 6J. Substituting into AMD 
produces the bound 

m ( m _l)ffll|| F ( )|| ||F(0)||^ (35) 
on the trace norm of the last term in (jHEJ)- 



5 Accuracy of the TDHF approximation 

We have shown that if D is a density operator on Fe of the form D = @D n such that 
J2n 2 Ti(D n ) < oo, then 

U ( t 1 \[V h2 M2{s)]. 1 ds , 

where Mi and M2 are the one-particle and two-particle reduced number density operators for a 
system which evolved under the dynamics ((TJ from the initial state D. In this section we will 
compare Mi(t) to the solution of the TDHF equation 

F(t) = U^FiO) - j:£ UP.IV, F(sr 2 2A 2 ] :1 ds 

F(0) = Mi(0) (36) 

whose initial condition is jVi(O). When the initial state D is a Gibbs density for noninteracting 
fermions, then the distance between Mi(t) and F(t) in trace norm can be controlled at short 
times fTheorem l5.1|) . In the mean field limit, F(t) is an asymptotically accurate approximation 
to Mi(t) at all times t, provided that the initial states are Gibbs states for noninteracting 
fermions (Theorem 15. 2|) . 

We state and discuss Theorems 15.11 and 15. 21 before going on to prove them: 

Theorem 5.1 Let G be the density operator on Fh of a Gibbs equilibrium state for noninter- 
acting fermions, as in Proposition ^. 6 A LetMiit) andM2(t) denote MiiWtiG)) and M2{Wt(G)) , 
respectively, where Wt is the dynamics with two-particle interactions defined in \2ty) . 

Let F(t) be the solution of the TDHF equation 1136)) . 

Let t denote (2||V||||.A/i||i)- 1 7i. Then 

\m)-F(t)i < l(vh) 2 (37) 

for t < r. 
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This theorem implies, for instance, that 



- F(r/2)l 



3 

< -. 

~ 2 



This is remarkable because there are about |||A/i||f interactions driving the dynamics and 
the error could be much larger prima facie: it could be proportional to ||A/i||i. Unfortunately, 
the bound (|37|) on the error of the TDHF approximation is valid only when t < r, which is 
inversely proportional to ||.A/i||i, and we have no explicit bounds for larger t. If ||A/i||i is too 
large, the bound (|37j) is useless, for then it is valid for too short a time. 

The time-of- validity of (J37j) ends up being inversely proportional to ||A/"i||i because the number 
of two-particle interactions is proportional to ||A/i||f, and none of these interactions is weaker 
than any other a priori. In the thermodynamic limit, where ||A/i||i — > oo with constant spatial 
density, the total interaction energy grows like ||A/i||i (rather than the square of ||JVi||i) if the 
interaction potential is short-ranged. We do not know how to derive the TDHF equation in 
the thermodynamic limit, but we can derive it in the mean field limit, where the strength of 
the interaction is scaled in inverse proportion to ||A/i||i. 

For each value of the parameter A > 0, consider the Hamiltonian 

H x = ^2(j\L\i)a]ai + \ ^{kl\V\i])a\a\a j a i . (38) 

If the initial state is given by a Gibbs equilibrium density D x for noninteracting fermions, then 
Proposition 13.11 implies that all reduced number density operators exist and satisfy 



l<i<j<m 

pt m 

J A 1 



AC(0) = Af m (D x (0)). (39) 

We will show that if A is inversely proportional to ||Ai ||i, then A/"„(t) is close to F x (t)® m m\A m 
in trace norm, where 

/•t 

F x (t) = U^F x (0) - % - l4? s [V, F x (s) m 2A 2 ]. 1 ds 

"Jo 

F A (0) = A/1 A (0). (40) 

Theorem 5.2 Let |Z) A } A>0 be a family of Gibbs equilibrium densities for noninteracting fermions, 
as in Proposition ^. <X with 

limsupA||A/'i( J D A )||i < oo. 
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Let M^{t) be the solution of let F x (t) be the solution of the TDHF equation ^Oj ) and let 
F*(t) = F x {t)® m m\A m . Then 



Km\\M^t)-^ m (t)\\J\\M?\\? = 

A — >() 

for allt > and all m e N. 

We will derive Theorems 15 . 1 1 and 15 . 21 from Lemma f5. II below. 

Let M m {t) be as in Proposition 13.11 and let T m {t) satisfy the TDHF hierarchy. In the 
hypotheses Theorems 15.11 and 15.21 we suppose that -F(O) = J\f\(D(0)), but for now let us only 
assume that 

||F(0)||<1 and 11^(0)11! = ||MP(0))||i. (41) 

The trace norm of Af\(t) is independent of t, and we shall denote it simply by Assuming 
flUJ, t ne bound (|33j) is itself bounded by 

m ( m _l)2M \\AfArt. (42) 
n 11 111 

Subtracting equations (jHSJ) from equations (fHU|) and using (|4*2*j) leads to the estimates 

\\Kn(t) ~ Fmit)^ < \\Kn(0)- ^(0)^ 

m(m — 1) 



+ m[m h - \\v\\(2Wi\\T + Wmh)t 

/ ||A/ m+ i(s) - ^^(s)!! ds 
Jo 



2\\V U rt 

+ m 



h 

Iterating this estimate n times, one obtains 

IKW-^WIL < ±a m J m+j r 1 )cv + ± b ^( m+j - 1 )cH'- 

j=o \ 3 / j =0 J + V 3 J 

(m + n)\ f l f h f tn 

+ Cn ~( TTT / / "/ I Wm+n+l (s) - ^m+n+l(s) || ^sdtn ■■ ■ dt x 

(43) 

with C = 2\\V\\/h and 



a. 



11^(0) -^(0)^ 



b m = m{m H 1) ||y||(2||AM.||;" + ||ACII.i.)- «in 
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To make use of these estimates we need some control over the size of the integrand in (J43|) . We 
will assume that there exists a constant B such that 

IIAUi < ^IIMIir (45) 
for all m. Then the last term on the right hand side of is bounded by 

'm + n — I s 



m — 



1 ^(s + ijllMlpVTV' 



which tends to as n tends to infinity if m fixed and C||jV"i||it < 1. Furthermore, assuming 
(Hoi , we can bound b m of flSJ) by C(B/2 + l)||A/"i||™m(m - 1) and we find the following: 

Lemma 5.1 Suppose that D is a density operator on Fe of the form D = ®D n , such that \4°V 
holds for all m e N. Let W t be as defined in H|) and let N m (t) denote J\f m (W t (D)). Let F(t) 
be the solution of a TDHF equation < f37)) whose initial condition F(0) satisfies fiTty , and let 
Tmit) be as in $3%). Then, with C = 2\\V\\/h, 



\\N m (t) - . 


T m {t)\ 


1 




|M| 


| m 
1 



+ iiMii^|;(».+i-i)g:j)mii.*r 



when C||jVi||i< < 1- 



Proof of Theorem 15.11 When G is the density operator of a Gibbs equilibrium state and 
F(0) =7Vi(G), then ||F(0)|| = \\Afi(G)\\ < 1 by Proposition O ||W ro (0) - F m (0) || t = for all 
m by Proposition 12.21 and B = 1 in (|45|) by Proposition 12.41 Upon simplifying the inequality 
in Lemma f5. 11 one obtains Theorem 15. II □ 

Proof of Theorem O For Af*(t) and F*(t) as in the hypothesis of Theorem 15.21 observe 
that 

||F(t)||<l and \\F{t)\\i = \WiW))\\i 
at all times t, and J\f^ satisfies f!45|) with B = 1. Thus, we may apply Lemma f5.1[ with a few 
changes: t and t + At may be substituted for and t, and C should be replaced by CX. These 
substitutions yield 

iP^Hr - g i^r^ I m-i J^ii^n^) 



(46) 
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for At < (CX\\Afi\\i) 1 . Since u = limsupCA||A/' 1 A ||i is finite by hypothesis, taking the limsup 
of both sides of implies that 

Um W^-^yil^X = V m £ N (47) 

holds at time s — t + At if it holds at time s — t and At < 1/u. Since (J47|) holds at s = 0, an 
inductive argument proves that it holds at all times s > 0. □ 

6 Appendix: the proofs of Propositions 12. 1L 12. 2L and 12.3 

6.1 Proof of Proposition 12.11 

We begin by proving 

Proposition 6.1 If D n is an n-particle fermionic density operator then 

\\D na \\ < 1/n. (48) 

Proof: Thanks to the convexity of the norm and the linearity of the partial trace, it 
suffices to prove ()48j) for fermionic pure states. Let EI denote the single-particle Hilbert space, 
and let \l/ be a unit vector in M^ n \ Since P^-\ is a compact Hermitian operator, there exists a 
unit vector u £ EI such that 

||P* :1 || = (u,Pu :1 (u)). (49) 

Let {4>j}j(zj be a basis of H containing u. For each subset s = {ji, . . . ,j n } of J, let ^(s) denote 
one of the two Slater determinants that may be formed from the vectors <\>j x , . . . , <fij n (the two 
choices differ only in sign). The set of vectors ^(s) is an ortho normal basis of El( n ) and so 
^ = ^ s (\E'(s), \E')\[ / (s). By definition of the partial trace, (u, -FW(it)) equals 

(« ® 4>h ® ■ ■ • ® &„_i, Mu ® (f>n ® • • • ® (50) 

ii,...j n -ieJ 
From and flSUJ) 

ji,...,j„-ieJ s 
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But 

1/n! if s = {u,<f> jl ,...,<f> jn _ 1 } 



whence 



otherwise 



iip* :1 ii < j:\(ns)^}\ 2 ^^ = -Ek*( s )>*>i 2 = 1 

z — ' n! n ^— ' n 

s s 

by dSIJ). □ 

Now we return to the proof of Proposition 12.11 

Let D be a density operator on Fe that commutes with the number operator N and satisfies 
Tt(ND) < oo. Then D has the form JUJ) and £>Tr(.D n ) < oo, so that Ni{D) = J2 nD n-i is 
defined and 



\\Ni(D)\\ < £n||£) n!l ||. 



n=l 



By Proposition 16. II the operator norm of D n is less than or equal to Tr(D n )/n, whence 

oo 

||M(-D)|| < ^TnTr(D n )/n = Tr(D) = 1. 



n=l 



6.2 Proof of Proposition 12.21 

From (|221), with p 0) as in pi . 



s:N(s)>n K W 7 

Substituting the expressions (fT5j) and (|T2|) for (P*( s )). and collecting terms, we find that 
A^.(G) = £ P(s)(AT(s)P* (s) ) : f n nL4 n 

s:N(s)>n 

= E [ E P ( S )] ® • • • ® P ^> !A n , (52) 

distinct S3jl • .ijn 
il,— JnSJ 

where iV(s) is the size of s. The sum in (I52J) is made over distinct ji, . . . ,j n since (P^ ® • • • ® 
P^.^)A n equals the zero operator if j r = j s for any r ^ s. In particular, 

= E[E p ( s )H = EfO")^- (53) 
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Thus, 

n 

Ni{G)® n n\A n = \{p{ji)(P<t> h ®'--®P^ n )n\A n . (54) 

distinct i— 1 

The sum in (p>f)) is again restricted to distinct ji, . . . ,j n because of the presence of the anti- 
symmetrizer A n . Substituting (J2Tj) into (jo^j) yields (Jo"2|) . proving the proposition. 

6.3 Proof of Proposition 12.31 

Suppose T is a Hermitian trace class operator. There exists an orthonormal basis {&j}jeJ of EI 
such that 

T = with = H T llr 

i j 

The operator T® n n\A n is diagonalizable with respect to the basis of Slater determinants formed 
from n distinct members of {e 3 }. Indeed, 



T® n n\A n (Vn). A n {e h ® ■ ■ ■ ® e jn )) = Vn\A n n\T^ n {e jl <g> • • • ® e jn ) 

n 

= (n\ J] A is ) Vri. A n ( eil ® • • • ® e in ). 

s=l 

The trace norm of T* x,n n!y4 n is the sum of the absolute values of its eigenvalues, whence 

n 

WT^nlAnl = Yl IIl^l- ( 55 ) 

distinct S = \ 

Note that the sum in (J55J) is over ordered sequences ji, j 2 , • • • , jn rather than subsets {ji, . . . , j n }. 
But 

n n n 

e niA*.i < e niA*.i = iiem = ™ 

distinct S = l jli — ~,jn£.J S=l S = l j£j 

ji,---,j n eJ 

proving 
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